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Finite Volume Effects in Weak Hadronic Decays * 

C.-J.D. Lin^, G. Martinelli'^, C.T. Sachrajda^,M. Testa*'. 

'^Dept. of Physics and Astronomy, Univ. of Southampton, 
Southampton, S017 IBJ, UK 

^Dip. di Fisica, Univ. di Roma "La Sapienza" and INFN, Sezione di Roma, 
Piazzale Aldo Moro 2, 1-00185 Rome, Italy 

In this talk we discuss finite-volume computations of two-body hadronic decays below the inelastic threshold 
(e.g. K — > TTTT decays). In particular we show how the relation between finite- volume matrix elements and physical 
amplitudes, recently derived by Lellouch and Liischer, can be extended to all elastic states under the inelastic 
threshold. We also provide a derivation of the Liischer quantization condition directly in quantum field theory. 



1. Introduction 

Lattice QCD offers a natural opportunity to 
compute the matrix elements for K tttt de- 
cays from first principles. The main difficulties 
are related to the continuum limit of the regu- 
larized theory (the ultra-violet problem) and to 
the relation between matrix elements computed 
in a finite Euclidean space-time volume and the 
corresponding physical amplitudes (the infrared 
problem). The ultra-violet problem, which deals 
with the construction of finite matrix elements of 
renormalized operators from the bare lattice ones, 
has been addressed in a series of papers and 
we will not consider it further. The infrared prob- 
lem arises from two sources: 

• the unavoidable continuation of the theory 
to Euclidean space-time and 

• the use of a finite volume in numerical sim- 
ulations. 

One of the main obstacles in the extraction 
of physical amplitudes from lattice simulations 
stems from the rescattering of final state parti- 
cles in Euclidean space. The formalization of this 
problem, in the infinite-volume case, was consid- 
ered in ref. An important step towards the 
solution of the infrared problem has recently been 



achieved by Lellouch and Liischer Q (LL), who 
derived a relation, given below in eq.(|l^), be- 
tween the K —^ TTTT matrix elements in a finite 
volume and the physical kaon-decay amplitudes. 
For technical reasons the LL formula has been 
derived for a finite, fixed number of pion states 
under the inelastic threshold and for matrix el- 
ements at zero four-momentum transfer. In the 
following sections a different approach to the LL 
formula will be discussed which extends it to 
all elastic states under the inelastic threshold and 
momentum transfers different from zero. 

2. Physics in a Finite Cubic Volume 

In this section we will recall some aspects of 
the quantization in a cubic box, relevant in the 
discussion of the LL formula In partic- 

ular we will be interested in the structure of 
the zero-momentum, finite volume, energy eigen- 
states, Itttt, n)y, which can be excited from the 
vacuum by a scalar operator cr(x) i.e. for which 



V (7r7r,n|cr(0) |0) ^ 



(1) 



When analyzed from the point of view of angular 
momentum, energy-eigenstates in a cubic box are 
a complicated superposition 
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We are interested in the structure of states which 
are cubically invariant and contain an s-wave 
component (and hence satisfy eq.(|])). 

In the framework of quantum mechanics, finite- 
volume quantization formulas become exact for 
finite-range potentials and in the presence of an 
angular momentum cut-off, which, in its simplest 
form, assumes that only s- waves interact. Un- 
der these conditions, the allowed values of the 
radial relative momentum /c of a two particle 
state, related to the center of mass energy E as 
E ~ 2^m^ + fc2, are quantized as follows 0: 



1. fc obeys the equation (see sec.^) 



(3) 



where n is a non-negative integer, (5(fc) is the 
infinite volume s-wave phase-shift, q = 
and 



tan <^{ci) 



(4) 



^oo(i;92) 

Zoo(s;g2)^^ ^ {n^ - q^Y' .(5) 



2. — — (free spectrum) if at 

least two p and p' exist with p' ^ = not 

—71 i^n — n — n 

related by a cubic transformation. These 
states are non physical and their existence 
is a consequence of the angular momentum 
cut-off. 

States of type ^ have a non-zero s-wave compo- 
nent (r), undistorted by the presence 
of the boundary compared to the infinite volume 
s-wave function x^/'^"-""^^' (^r). Therefore inside 
the volume V, 



^c{En) 



(6) 



States of type g cannot be simply plane- waves of 
the form e*-'> because of the presence of the 
interaction. However the combination 4>{x) — 
is a solution of the Schroedinger 



equation even in the presence of the potential, be- 
cause (/)(a:) does not contain an s-wave component 
and all other angular momenta are not interact- 
ing The spurious states (j){x) have a non-zero 
cubically-invariant projection, but do not project 
on s-wave. 

In conclusion, locality and the scalar character 
of a{x) imply that 



(0| cr(0) Itttt, n; Spurious) y = 



|a(0)|7r7r,n)^= / ^^^^(r)*^ Jx) 

JV 



(7) 



d^xS'(r)^'^=-™""=(r) 



(8) 



where S{r) is the coordinate representation of 
a{x). S{r) is a function of r only, so that it se- 
lects the s-wave component x^/^^-""'^' (?■) in the an- 
gular momentum expansion, eq.(H). S{r) is also 
localized inside V, thus justifying the last step in 
eq.(||). With the definition 



a{En) = {0\a{0) Itttt, En) 
d3a;^(r)*^=-™""=(r) 



we therefore have 

\aiEn)\^ ^ c{E.n)ma{0)\7T7T,n)yf 



(9) 



(10) 



After these preliminaries we are ready to discuss 
the LL proposal, which consists in tuning the vol- 
ume V so that one of the first seven excited two- 
pion state is degenerate in energy with the kaon 
state (for n = l, L^ib-^Q Fm) and then using 
the LL relation which connects finite and infinite 
volume matrix elements 

\{7:7:,E^mK\nw{0)\K)\^ = (11) 
= V' \v {tttt, E\ HwiO) \K)y\' X 
x8TT[q(j>'{q) + k6'{k)] . 



In eq.(|ll|) \TnT,E)y denotes a finite volume two 
pion state with zero total momentum and "angu- 
lar momentum" , whose energy E is to be chosen 
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among the solutions of eq.(g), while \K)y denotes 
a finite volume kaon state with zero momentum. 
Both states are normalized to 1. The LL formula 
has been derived |^ for a large enough volume, 
n = 1 ^ 7 and A£; = AP = 0. 

Eq.jl^) shows that c{En), defined in eq.(||), is 
precisely the LL proportionality factor. 

3. The Nature of the LL Relation 

In order to relate the states at finite and infinite 
volume we consider the correlator 



(fx {a{x,t)a{0))y 
(27r)3 f dp^ dp^ 



(12) 



x|(0|a(0)b^,p^)P = 



-(ui+W2)t 



2(27r)3 
dp^ dp^ 



dEe 



|ct(0) \TnT,E)\^ X 



= 2(^ I ^e-^*l(0k(0)k7r,i?)|^fc(£;). 
On the other hand we could proceed differently 
d^x{a{x,t)a{0)) = (13) 



y^|(0|a(0)|7r7r,i 



|2 



V / dEpv{E)\{0\a{0)\TnT,E)y\\ 



where only states containing an s-wave compo- 
nent contribute. In eq.(p^ Pv{E) denotes a func- 
tion to be determined, which provides the corre- 
spondence between finite volume sums and infi- 
nite volume integrals. In many cases, for example 
in one dimension, pv{E) can be identified as the 
density of states at energy E. In section ^ we 
show that, also in the presence of cubic boundary 
conditions, pv{E) is given by 



dn q<P'{q) + k5' [k) ^ 
= dE= ^ 



(14) 



with exponential precision in the volume. The ex- 
pression in eq. (|lj) is the one one would heuristi- 
cally derive from a nai've interpretation of pv {E) 
as the density of states, as seen from eq. (||). 

Comparing eqs.dl^) and (^ we get the corre- 
spondence 



Itttt, E) -ifi^ 47ri 



IVEpvjE) 
HE) 



Tnr,E)^ 



In a similar way one gets 



\K,p^O) V2mV\K,p = 0) 



(15) 



(16) 



from which it is easy to recover the LL rela- 
tion, eq.(pT|), without any restriction on the four- 
momentum transfer. 

Although the present approach appears super- 
ficially to be equivalent to the one of ref. , there 
is an important difference in the two derivations. 
The result of ref. § was obtained at a fixed value 
of n and therefore at a fixed volume V, tuned 
so that rriK = En, with n < 8. We, on the 
other hand, derived the same result at fixed en- 
ergy E, for asymptotically large volumes V. This 
implies that, as 1^ — s- cx), we must simultaneously 
allow n — !■ oo. A question which arises naturally 
at this point is what is the relation between the 
two approaches? The answer requires a more de- 
tailed discussion, developed in the following sec- 
tion, where it will be shown that the constraints 
of locality allow us to establish eq. (|lj) with ex- 
ponential accuracy for elastic states under the in- 
elastic threshold. 

4. Summation Theorems, Locality and the 
LL Formula 

Locality has already been an important ingre- 
dient in establishing eq.(|o|). In this section we 
will discuss another important and more subtle 
role of locality. Our approach to the LL formula, 
outlined in the previous section, is based on the 
identification 



E I (01^(0) I 



dE\a{E)\\''^* 



(17) 



4 



of finite and infinite volume correlators. We are 
therefore naturally led to the question of how well 
a sum may approximate an integral. We will show 
that the key ingredient to answer this question is 
again locality. 

We start with the example of a simple summa- 
tion theorem. Take any f{x) and g{x) exponen- 
tially decreasing at large |^| and compute their 
Fourier transforms 

f{p) ^ Jdx fix)e't-^ , g{p) = J dx g{x)e'f^. (18) 

V V 

OO CO 

In eq.(p^ the integrals can be taken on a finite 
volume V or over all the three-dimensional space, 
the difference being exponentially small in V. If 
we define 



2tt 

p = — n 
L - 

we have 



(19) 



Kx)g*(x)dx -■ 



(27r)3 



f{p)g*{p)dp 



= mr{2L)dx ^ 1 ^ fiPj9*iPj (20) 



'V 

and therefore 
1 



(27r)3 



f{p)g*{p)dp 



V ^ 



/(« )g*(« ), (21) 



again with exponential accuracy in V . If the sup- 
port of /(x) and g(x) is entirely contained inside 
Y , eq.(pi|) becomes exact. 

Relations similar in nature to eq.(^l]) can be 
obtained in quantum mechanics and field theory. 
In fact, since (y(x) is a local operator, (t(0) |0) is 
a localized state which does not differ much from 
the vacuum state away from 0. This is a con- 
sequence of clustering, which guarantees that, in 
the absence of massless particles, the probability 
of finding particles at a distance r away from the 
origin in the state (t(0) |0) decreases exponentially 
like e"^™"''. As a consequence, if V is greater that 
the localization volume, we can write, at t = 0, 



0)cr(0))y = |(0|cr(0)|7r7r,n) 



dE\a{E)'^ ^ d^x{a{x,0)a{0)) (22) 



with exponential accuracy in V. The second and 
third lines of eqs.(^) were obtained by inserting 
the complete set of energy eigenstates on a finite 
and infinite volume respectively. The result is the 
same apart from the exponentially small pertur- 
bation at the boundary. In other words, we can 
introduce states 



(23) 



where {x) are the finite volume center of mass 
wave functions. |7r7r,x)y represents a state with 
P = 0, in which the two pions are localized a 
distance x apart. If the distance x is inside V, 
the same state is exponentially well represented 
by 

|7r7r,x) = (27r)3/2 / di; Itttt, *f *(x) (24) 
Je 

where "^^{x) are the infinite- volume center of 
mass wave functions. 

By virtue of the cluster property mentioned 
above, we know that 

F(x) = (7r7r,x|CT(0) |0) (25) 
is a localized function of x, so that 
P(x) « y(7r7r,x|cr(0)|0) (26) 
with exponential precision and 



d^x{a{x,0)a{0))v 



0\a{0)\TnT,x)v\^dx-- 



(27) 



|(0|cr(0)|7r7r,a;)|^da; : 



dE\a{E)\'' 



d^x{a{x,0)a{0)) 



confirming cq. (p^) . Eqs.(p2|) and (||) express the 
important property that the matrix elements of a 
local operator are smooth in energy, in the sense 
that the sum over discrete energy levels repro- 
duces, with exponential accuracy, the correspond- 
ing integral 



POO 

/ dE\a{E){'^VJ2- 



c{E^) 



(28) 
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Introduction of time dependence does not sub- 
stantially change this result. We have 



dPx{a{x,t)c7{S)))v = 



(29) 



\2-E„t 



= y^|(0|(T(0)|7r7r,n)ype- 

1 1 cr (0) I TTTT, x) e (tttt, X I fj (0) 1 0) = 
= / F*{x){e-~"'F){x)dx^ [ F* {x)F{x,t)dx, 



where the two-body hamiltonian H is related to 
the second quantized hamiltonian, H, through 
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= \/V^^|7r7r,n)ye"-^**|jJx) . (30) 



F{x,t) as defined in cq.(p9[) has a time evolution 
characteristic of a diffusion process with hamil- 
tonian H. Therefore both F{x, t) and the corre- 
lator Jy (Px{a{x,t)a{0))v are insensitive to the 
presence of the boundaries, up to exponentially 
small terms, for a long time (t ~ m.,^L), and we 
have 



dE\a{E)Ye 



\a{E„ 



-E„t 



c{E„) 



(31) 



with exponential precision in V . 

The basic tool for relating integrals to sums is 
provided by the Poisson identity 



+ 00 



(32) 



1 — — 00 



Eq.(32), together with the substitution x — > 
h{E,L), where h{E,L) is defined in eq.(||), gives 



E 



a{E)\'e~ 


Et _ 


|fT(S„)|'e 


-_E„t 


dh(E,L) 




BE 


Er^ 



(33) 



where 



f-OO 

Q{L,t) = Y^ / di;|cr(£;)|2e-^*e2*'*('=)e2*''^(«).(34) 



Consider the large L behavior of Q(L,t). The 
techniques of asymptotic analysis [||, together 
with the fact that 0(g) — > oo when L — s- oo at 
fixed k, suggest that the large L behavior of each 
term of the sum over / in eq.(Q) is dominated 
by the critical points of ^(g), i.e. the points 
at which (f)[q) has a vanishing derivative or the 
points where and cr{E) are not continuous 
or differentiable. 5 = is a critical point of (l){q) 
as follows from 



2„3 



g«0 



(35) 



so that, in absence of further critical points we 
expect that 



Q{LA) 







(36) 



as a power of j- depending on the threshold be- 
havior of <7{E) which, for local observables, may 
be make arbitrarily small |^ . 

A comparison between eqs.(31) and ( ^ ) then 
leads to the identification 



c{En) ^ V 



dh{E,L) 



dE 



(37) 



and gives the LL relation with exponential accu- 
racy, for energies between the elastic and inelastic 
thresholds. 

Although a rigorous mathematical proof that 
g = is the only critical point of (/)(g) in the com- 
plex g-plane is lacking, we can provide a simple, 
indirect argument that this is likely to be the case. 
For this we start applying eq.(^ to the case of 
two functions f{p) and g{p) which only depend 
on p = \p\ 

1 1"^ 1 

^ f{p)g*{p)p'dp = - ^ f{Pn)9*{Pn) = 



{P } 



= E '^nfiPn)9*iPn) , 



(38) 



where Vn is the number of integer vectors with 
given The last sum in eq.(|3^) is performed 
over the distinct values of p„. On the other hand 
the PnS are the solutions of 
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because, from eq.(^ we have that 
tan (t){q) — ^ q — \n\ 



(40) 



where n is any vector with integer components. 
We can now use the Poisson Identity eq. ( |3^ ) with 
X — > (/)((j)/7r which, together with 

47r2 

<l>'{\n\)\ = —n' (41) 

and after multiphcation by p'^ /2n^f{p)g*{p) and 
integration over p, gives 

3^/(0)5*(0) + -^ J2 ^nf{Pn)9*{Pn) = 



^/(0)5*(0) + ^E^"/(?'")ff*(^^») 



(27r)3 

y- 



V 

{p>.} 
f{p)9*{p)dp 



p'f{p)9*{p)e'''^^^^dp. (42) 



Eq.(^), together with eq.js^) shows that 

1 poo 

EA/ P'm9*{p)e'''''^^^dp (43) 



is indeed exponentially small with power correc- 
tions concentrated at threshold. This is in com- 
plete agreement with the heuristic analysis pre- 
sented after eq. ( p^ ) . This argument strongly sug- 
gests that the only critical point of (/)(g) is at 
q = 0. In fact if (j){q) possessed other critical 
points in the complex g-plane, these would show 
up in eq.(^), through a large L behavior different 
from the one just found. 

5. Finite Volume Quantization in Quan- 
tum Field Theory and the Effects of In- 
elasticity 

The approach presented in previous sections is 
based on the property that local correlators may 
be expressed, with exponential accuracy, both as 
a sum or as an integral over intermediate states. 
This may appear to be a great difference of our 
approach compared to that of ref. We will 
argue, however, that the validity of eq.(H) also re- 
quires the volume to be sufficiently large for the 



Fourier series to be equal to the infinite-volume 
energy integrals, up to exponential corrections. 
Since the derivation of the LL relation in ref. ^ 
relies on this quantization formula, we conclude 
that the conditions on the volume for the appli- 
cability of this relation are equivalent in the two 
approaches. 

The validity of eq.(||) in Relativistic Quan- 
tum Field Theory (RQFT) has been discussed in 
ref. 0. In this section we describe a different 
approach which helps to clarify the size of the 
corrections to eq.(^ due to the presence of an 
inelastic threshold, Ei„ 

The concept of a wave function in RQFT is an 
approximate one. The object closest to a wave 
function is the Bethe-Salpeter (BS) wave func- 
tion. In infinite-volume, for an incoming state 
with total momentum zero, the t — BS wave 
function is defined as 



$fc(ai) = (O|(/.(x,O)0(O)|fc,-fc),„ = 
= E(O|0(2,O)|n)(n|0(O)|fc, -fc), 



(44) 



where 4'{x) is an appropriately normalized pion 
field. In eq.(^) we can separate the contribution 
of single pion states, vE'fe(x) and that of multipion 
states, Ik{x), as 



$fc (x) = *fc(x) -I- Ik{x) 
where 

d^p 



(45) 



(p|<^(0)|fc,-fc)„e%^ (46) 



{3^} 

+ ••• (47) 

where the second sum runs over three pion states 
and the ellipses represent the contribution from 
states with a higher number of particles. As will 
be shown later, Ikix) vanishes exponentially with 
X and will be neglected for the moment. We can 
further separate, in ^'fc(a;), the connected part 



*fc(^) 



(48) 



d^p 



{2TTf2Ep 



{p\ m \k,^k)i 
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and parametrize (p|0(O)|fc, — fc)™"' 
(p|<^(0)|A:,-fc)r" ^ ^ 



4Ek Ep-Ek- le 



(49) 



The off-shell scattering amplitude A4 becomes the 
physical one 0, M{k^p), when p — k 



— i k 
We therefore have 

ik-x I 



(50) 



(51) 



(27r)32i;p4i;fc(£;p-^fc-i£) 



The projection of 'I'a.. (x) over the s-wave is 



\ s — wave 



(r) 



(52) 



(53) 



{2Tr)32Ep AEkip"^ - k^) pr ' 
where we used the identity 

— — = in6{x)+V- . 

X — ie X 



From eq.(p2|) we obtain the large r behaviour of 

^ k] s— wave (^) 



(,.) = J_ ( sinfcr + l(e*2'5(fe) - l)e^'=M 
fcr \ 2i y 



kjs—wavey' ) 

As discussed in section ^, the relative momenta k 
allowed in a box are determined by the condition 
that the s-wave projection of the finite volume BS 
wave function, ^'^\s-wave{'''), is not deformed by 
the presence of the boundary. <&^(a;) is given by 
an expression similar to eg. (^) , where the sum 
over intermediate states runs over finite volume 
energy eigenstates, so that, neglecting for the mo- 
ment the contribution from multipion states, we 
have 



^fe \s — wave{f) 



(55) 



1 (Ep^ + Ek)M sin p„r 



V ^^2Ep^ 4Ek{pl-k^) pr,r 



(54) 



where the sum is over single particle momenta 
given by eq.(|l^). The quantization condition re- 
quires therefore the identity of eqs.(52) and (^. 
We stress that in the finite volume expression, 
eq. (p5|) , disconnected terms and ie's are absent 
because the eigen-momenta k do not coincide 
with any of the p„. These terms will appear as V 
grows, as shown in a moment. 

In appendix C of ref. |^ we prove the summa- 
tion formula 



-y 

Vf^k'^ 



fi.pl) 



(56) 



(27r)3 



c/(fc') 



where 



^oo(l,g') 
2tt^L 



(57) 



As before we keep the interaction only for the s-wave. 



Eq.(|5^) is valid up to exponentially small cor- 
rections, provided the Fourier transform of /(p^) 
vanishes exponentially with r and, together with 
eqs.(|5|) and (p2|), gives at once the Liischer quan- 
tization condition, eq. (^) , with exponential preci- 
sion in V. The size of the corrections are de- 
termined by the precision with which the sum in 
eq.(p5|) is able to reproduce the corresponding in- 
tegral, a condition analogous to the one which was 
underlying our derivation of the LL formula. It 
appears therefore that the two formulas have the 
same conditions of applicability. In particular, a 
possible source of concern is the (practically rele- 
vant) situation in which the quantization volume 
allows only very few (perhaps two or three) elastic 
states under the inelastic threshold. In this case, 
even though the finite-volume effects are expo- 
nentially small, one is working in a fixed volume 
and in order to establish that the corrections are 
indeed negligible one needs an estimate of such ef- 
fects. These corrections are small if '^k(j')\s-wave 
is well approximated by its asymptotic expres- 
sion, eq.(|5^), at r w L. There are essentially 
two types of exponentially small corrections. The 
first type is analogous to the corrections which 
would be present, even in quantum mechanics in 
the presence of an exponentially decreasing po- 
tential, rather than a finite range one. The sec- 
ond type of corrections, characteristic of RQFT, 
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are due to the existence of inelasticity and are re- 
sponsible for the failure of eq.(||) above Ein. They 
come from the term Ik{x) in eq.(^) and a cor- 
responding crossed contribution in '^k{x)- The 
matrix element (37r|0(O)|fc, —k)in in eq.(p7|) has a 
singularity when the virtual pion described by (f) 
reaches the mass-shell. This singularity is com- 
plex for 2Ek < Ein , implying an exponential van- 
ishing of Ik{x)- For values of 2Ek closer to E^ 
the singularity is closer to the real axis, and the 
range of the exponential increases. Therefore the 
quantization condition is, in general, affected by 
the presence of inelasticity. In the case of pions, 
however, there are indications that inelasticity is 
negligible up to around 1 Gev . This shifts the 
singularity further from the real axis, making the 
finite-volume effects less important and allowing 
the LL strategy to be applicable to — > tttt de- 
cays without sizeable corrections. 

6. Conclusions 

We have reviewed finite- volume effects in K ^ 
TTTT decays, examining the conditions under which 
the LL formula is valid. We were able to re- 
move the restrictions of zero momentum transfer 
and n < 8, present in the original LL deriva- 
tion, allowing the infinite volume limit to be 
taken. We have also presented an alternative sim- 
ple derivation of the Liischer quantization condi- 
tion in Quantum Field Theory, which allows us to 
examine the general consequences of inelasticity 
on the quantization formula. For the particular 
case of — » TTTT decays we conclude that these 
corrections are not important. 



M. Bochicchio et al., Nud. Phys. B262 (1988) 
B262; 

L. Maiani, G. Martinelli, G.C. Rossi and 
M. Testa, Phys. Lett. 176B (1986) 445 and 
Nud. Phys. B289 (1987) 505. 
C. Bernard, T. Draper, G. Hockney and 
A. Soni, Nud. Phys. B.Proc.Suppl. 4 
(1988) 483. 

C. Dawson et al., Nud. Phys. B514 (1998) 
313. 

L. Maiani and M. Testa, Phys. Lett. B245 
(1990) 585. 

L. Lellouch and M. Liischer, Commun. Math. 
Phys. 219 (2001) 31. 

C.-J.D. L in, G. MartineUi , C.T. Sachrajda, 
M. Testa, |hep-lat/0104006| , Nud. Phys.B, to 
appear. 

M. Liischer, Commun. Math. Phys. 104 
(1986) 177; Commun. Math. Phys. 105 
(1986) 153; Nud. Phys. B354 (1991) 531; 
Nud. Phys. B364 (1991) 237. 
N. Bleistein, R.A. Handelsman Asymp- 
totic expansions of integrals (Dover Editions, 
1986). 

See for example D. Morgan and M.R. Pen- 
nington, in The second Daphne Physics 
Handbook, edited by L Maiani, G. Pancheri 
and N. Paver, INFN-LNF. 



Acknowledgements 

M. Testa thanks the organizers of LHP2001 for 
the hospitality in Cairns. D. Lin, C. Sachrajda 
and M. Testa thank the Institute for Nuclear The- 
ory at the University of Washington for its hospi- 
tality and the Department of Energy for partial 
support during the completion of this work. This 
work was supported by European Union grant 
HTRN-CT-2000-00145. 



REFERENCES 



